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Abstract. In this paper we study the ground state phase diagram of a one-dimensional t — U — J model, 
at half-filling. In the large-bandwidth limit and for ferromagnetic exchange with easy-plane anisotropy, 
a phase with gapless charge and massive spin excitations, characterized by the coexistence of triplet su- 
perconducting and spin density wave instabilities is realized in the ground state. With reduction of the 
bandwidth, a transition into an insulating phase showing properties of the spin-^ XY model takes place. 
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1 Introduction 

Superconductivity near a magnetic instability is a topic 
of increased current interest in condensed matter physics. 
Magnetically mediated Cooper pairing near the antifer- 
romagnctic instability is widely discussed in the context 
of superconductivity in copper-oxide systems More- 
over, the discovery of Triplet Superconductivity (TS) in 
Sr2RuOi 121 and the recent discovery of coexistence of 
the TS phase with ferromagnetism in UGe2 PI , U RhGe ^ 
and ZrZni [S] has triggered an increased activity in stud- 
ies of correlated electron models showing close proximity 
of triplet superconducting and ferromagnetically ordered 
phases |tilVl8l9l1()l11ll2l1,Sini . 

Another group of unconventional superconductors with 
close proximity of magnetic and superconducting ordering 
belongs to the {TMTSF)2X family of quasi-one-dimen- 
sional conductors (Bechgaard salts) |15ll6j . Growing ex- 
■ perimental evidence has been collected in the last few 
; years, indicating that {TMTSF)2ClOi and {TMTSF)2PFq 
, under pressure are triplet superconductors |17j . Most in- 
' teresting is the phase diagram of {TMTSF)2PFq which 
' shows a spin-Peierls (SP) phase in the ground state at 
atmospheric pressure. Increasing pressure leads first to a 
transition from the SP phase into a spin density wave 
(SDW) phase, and finally to the suppression of the SDW 
ground state in favor of superconductivity ^Sj. Recent 
detailed experimental studies of the phase diagram of the 
(TMTSF)2PFq compound indicate the possibility of a 
coexistence regime between SDW and superconductivity 



Send offprint requests to: A. Schadschneider (as@thp.uni- 
koeln.de) 

^ Permanent address: Andronikashvili Institute of Physics, 
Georgian Academy, Tamarashvili 6, Tbilisi 380077, Georgia. 
Electronic address: japa@iphac.ge 



|19| . Although the very presence of a homogeneous coex- 
istence phase in the phase diagram was questioned in the 
more recent publication [201, the SDW-SC competition is 
common in organic materials |21| and therefore models 
of correlated electrons exhibiting such phases are of great 
interest. 

Various models of strongly correlated electrons show- 
ing close proximity of magnetic (ferromagnetic) and su- 
perconducting (triplet superconducting) phases have been 
subject of intensive research in attempt to construct a the- 
oretical model for new superconducting materials. Usually 
these models are based on some extensions of the Hub- 
bard model. In particular, several extended versions of 
the repulsive Hubbard model have been employed for a 
long time as standard models for metal-insulator transi- 
tions, antiferromagnetism and high- Tc superconductivity 
1^ . At the same time, the Hubbard model in the case of 
sufficiently narrow band and/or low doping is a standard 
model for metallic ferromagnetism of itinerant electrons 

m 

Taking into account the experimentally observed easy- 
plane anisotropy of the spin exchange |24| in some of 
these materials, Japaridze and Miiller-Hartmann |25| pro- 
posed a rather simple extension of the Hubbard model 
with transverse {XY-tyj>e) anisotropy as a suitable ap- 
proach to such systems with coexisting orders. Indeed, 
this model was shown to exhibit an extremely rich weak- 
coupling phase diagram. In particular in the case of a half- 
filled band the weak-coupling ground state phase diagram 
consists of two insulating antiferromagnetic phases with 
easy-plane anisotropy and a spin gapful metallic phase 
with an identical decay of the triplet superconducting and 
spin density wave (SDW'^-') instabilities. Strong evidence 
for the presence of an additional transition into a ferro- 
magnetic XY phase has also been given |25| . 
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The model describes a system of itinerant electrons 
with transverse spin-exchange interaction between elec- 
trons on nearest-neighbor sites. The one-dimensional ver- 
sion of the Hamiltonian reads: 

+ uY,p^{n)p,in) + hLj2^s+S~^, + h.c.) . (1) 

n n 

Here cjj (c„q,) is the creation (annihilation) operator 
for an electron at site n with spin a, p,a{n) 



cr^g'^n where tr* {i = x,y,z) are the Pauli 



matrices. Below we restrict our consideration to the case 
of repulsive on-site interaction [/ > while the sign of the 
exchange interaction is arbitrary. 

One can easily verify that besides the obvious U{1) 
spin symmetry in the half-filled case the model is char- 
acterized by a SU{2) charge symmetry. An electron-hole 
transformation for one spin component interchanges the 
charge and spin degrees of freedom, and maps to the 
attractive Hubbard model with pair-hopping interaction 
|W7| . 

That the TS phase can be realized in ID correlated 
electron systems is well known from standard "g-ologi/^ 
studies The extended (U-V) Hubbard model with 
nearest-neighbor attraction {V < 0) has been intensively 
studied to explain the competition between SDW and su- 
perconducting instabilities in TMTSF compounds |29) . 
However, due to spin rotational invariance, in the extended 
Hubbard model the TS phase is realized only in the Lut- 
tinger hquid phase for \U\ < -2V \2mmi\ . where both 
charge and spin excitations are gapless. Singlet supercon- 
ducting (SS) and TS correlations show identical power-law 
decay at large distances and the TS instability dominates 
only due to weak logarithmic corrections [SJ • On the other 
hand, in the spin gapped phase U < 2V, the dynamical 
generation of a spin gap leads to the complete suppression 
of the TS and SDW instabilities. In marked contrast, the 
ferromagnetic transverse exchange between electrons on 
neighboring sites provides the possibility for realization of 
the SDW-TS phase in the case of gapped spin excitation 
spectrum |2S]. 

In this paper we study the model ^ in the case of a 
half-filled band using the DMRG techniques. We investi- 
gate the excitation spectrum of the system as well as the 
behavior of various correlation functions. Our numerical 
results confirm the predictions of a weak-coupling anal- 
ysis. In addition, we study in detail the strong-coupling 
sector of the phase diagram, focusing our attention on the 
ferromagnetic transition. With increasing transverse ferro- 
magnetic exchange this will reveal the possibility for two 
different scenarios of transition into the easy-plane XY 
ferromagnetic phase: in the case of weak on-site repulsion 
([/ < Uc — 2t) the first transition at J^y ~ ji'^y^ takes 
place from a spin gapped metallic phase into the insulat- 
ing SDW'^' phase with long-range order (LRO) and the 
latter becomes unstable towards the spin-flop transition 

into the ferromagnetic XY phase at Jxy ■ On the other 



hand, in the case of strong on-site repulsion, the metal- 
lic phase is absent and only the spin-flop transition from 
the SDW*-^^ phase into the ferromagnetic XY phase takes 

(c2) 

place at Jxy ■ The critical value of the exchange strongly 

fell 

depends on the on-site repulsion. At C/ = 0, Jxy = —it 



-At. In the case of weak U t these param- 



and Jxy 

eters are renormalized linearly in U, Jxy''' = at i7 ~ 

while for U ^ t we have jiy^"* ^ l/U. 

The paper is organized as follows: in the next section 
the weak-coupling continuum-limit version of the model 
is investigated. In Sect. |3I results of DMRG studies for 
chains up to L = 120 sites are presented. Finally, Sect. 01 
is devoted to a discussion and concluding remarks. 



2 The continuum-limit theory 

In this section we consider the low-energy effective field 
theory of the initial lattice model. Although this procedure 
has a long history and is reviewed in many places |32| , for 
clarity we briefly sketch the most important points and 
fix our notation and conventions. Considering the weak- 
coupling limit, |?7|,|J^j,| <^ t we linearize the spectrum 
and pass to the continuum limit by substituting 



z"i?„(a-) + (-z)"L„(a;) 



(2) 



X = nag, where oq is the lattice spacing, and Ra{x) and 
La{x) describe the Il{ight) and L(e/t) excitations with dis- 
persion relations E = ztvpp. These fields are assumed to 
be smooth on the scale of the lattice spacing and can be 
bosonized in a standard way |32j 



La{x) 



1 



V27rao 
1 



(x) 



— i\/i7r<pL 



V27rao 



(x) 



(3) 
(4) 



where ^R{L).a{x) are the Right (Left) moving Bose fields, 
carrying spin a. Next we define 



<t>R,a- 



(5) 



and introduce linear combination 

^c^{h + 4>i)l^, ^c^{B^ + ei)/V2, (6) 
= ((^T-'^i)/^. ^s = {e^~e^)/V2 (7) 

to describe the charge (c) and spin (s) degrees of free- 
dom, respectively. Then the Hamiltonian density of the 
bosonized model is given by 

H. = Tic + Tis 5 

■Hc^VcJ dx{^[{dx^cf + {dxdcf] 

M?. 



(v/8^^,)}, 

{\[{dx^sf + {dx^s? 

^cos(V8^V?s)|. 
an J 



Tis — J dx 



(8) 



(9) 
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Here we have defined 
2 (X," - 1) 

2(/C°-l) 
2ttM? 



9c 
9s 



ttvf 
1 



TTVp 



iU + Ji_), 
{U + Ji_), 



TTVF 
1 

TTVp 



{U + J±), 
{u - J±) , 



V, 



VF 



VF 



2t 1 



Jj_ 

2TTt 



(10) 

(11) 
(12) 

(13) 

(14) 



This mapping of the lattice electron system model onto 
the quantum theory of two independent quantum Bose 
fields described in terms of an "effective" sine-Gordon 
(SG) Hamiltonians © and Q will allow to extract the 
ground state properties of the initial model using the far- 
infrared properties of the quantum SG theory. 

The infrared behavior of the SG Hamiltonian is de- 
scribed by the corresponding pair of rcnormalization group 
(RG) equations for the effective coupling constants Kc{s) (0 
and Mc(5)(0 



dl 
dl 



= -2(i^,(,)(0-l)A-4(,)(0 



(15) 



where / = In(ao) and the bare values of the coupling con- 
stants are K,(,){1 = 0) = /CO^) and M,(,)(l = 0) = iV/O^^. 

The pair of RG equations (|15|l describes the Kosterlitz- 
Thouless transition pS] . 

The flow lines lie on the hyperbola 



4 (A' 



c(s) 



1 



is) 



Ks) 



/^'-4(i^,%-l)2-(M,%))2 

(16) 

and exhibit two different regimes depending on the rela- 
tion between the bare coupling constants |34j . 

Weak- coupling regime. For 2(/Cc(,s) — 1) > 
we are in the weak-coupling regime: the effective mass 
■Mc(s) — *■ 0. The low energy (large distance) behavior 
of the corresponding gapless mode is described by a free 
scalar field. 

The vacuum averages of exponentials of the correspond- 
ing fields show a power-law decay at large distances (77 = 
c,s) 



(17) 
,(18) 



and the only parameter controlling the infrared behavior 
in the gapless regime is the fixed-point value of the effec- 
tive coupling constants K*^^^^ = K^^s-^{1 = 00) determined 

from the Eq. (|16|l . Note that in the SU{2) symmetric case 
/i - and K*. = 1. 



Strong coupling regime. 



For 2{K'^ 



is) 



1) < 



is) 



the system scales to strong coupling: depending on the 
sign of the bare mass M^^^y the renormalized mass Aici 



is) 



is driven to ±00, signaling a crossover to one of two strong 
coupling regimes with a dynamical generation of a com- 
mensurability gap in the excitation spectrum. The flow o f 
l-^c(s) I to large values indicates that the Mc(s)C0sV8TTK(j) 
term in the sine-Gordon model dominates the long-distance 
properties of the system. Depending on the sign of the 
mass term, the field ip gets ordered with the expectation 
values 



s)> 



^7r/8A'c(s) 




(^^c%>0) 
iKs) < 0) 



(19) 



Using the initial values of the coupling constants, given 
in (|1()|) - (|13|1 . we obtain that flow trajectories in the charge 
sector (due to the S'C/(2)-charge symmetry) are along the 
separatrix gc = gu- Therefore, at 



U + Ja;y> 0. 



(20) 



there is a gap in the charge excitation spectrum [Ac 7^ 
0) and the charge fleld ipc is ordered with the vacuum 
expectation value 

(^c) - 0, (21) 

while a.t U + Jxy < the charge sector is gapless and the 
fixed-point value of the parameter K* is one. 

The U{1) symmetry of the spin channel ensures more 
alternatives. Depending on the relation between the bare 
coupling constants there are two different strong-coupling 
sectors in the spin channel. For 



U < min{0, Jxy} 



(22) 



the spin channel is massive {As 7^ 0) and the field (fs gets 
ordered with the vacuum expectation value 



while for 



iVs) - 0, 



Jxy < min{0, U} 



(23) 



(24) 



the spin channel is massive {Ag ^ 0), with vacuum expec- 
tation value 



{fs 



8R\ 



(25) 



In all other cases the excitation spectrum in the cor- 
responding channel is gapless. The low-energy behavior of 
the system is controlled by the fixed-point value of the 
Luttinger- liquid parameter K* = 1 + ^g*. In the par- 
ticular case of vanishing on-site interaction {U = 0) and 
antiferromagnetic exchange {Jxy > 0) one has to use a sec- 
ond order RG analysis to define accurately the fixed point 
value of the parameter Kg (for details, see Ref. |25p. 



4 C. Dziurzik et al.: Triplet superconductivity vs. easy-plane ferromagnetism in a ID itinerant electron system 



2.1 Order parameters 

To clarify the symmetry properties of the ground states of 
the system in different sectors we consider the following 
set of order parameters: 

1) the on-site density operator p{n) = p-\{n) + pi{n) 



p{n) 



Y^ici^aCn^a " 1) - \/^^: 



+ (-1)" sin( V2^(Pc) cos(y2^^,) , (26) 
2) the on-site spin-density 



In 



+ (-1)" cos(V27r/^c¥'c) sin (^2^^,) , (27) 
Sxin) 



2 / / 27f 

■ cos(v27rA^<(2s) sin ( \l ^-d. 



+ (-1)" cos(V27ri^e^,) sin ( \l^^s ) , (28) 




Sy{n) 



_j / 2tt 

COs{y^2TTKs(Ps) COS \ — t?, 

TTflo VV As 



+ (-1)" cos(\/27rA'cV'c) sm (^^j^^^ 



(29) 



and in addition we use superconducting order parameters 
corresponding to 

3a) the on-site singlet 



cos{^2tt Ks^Ps) exp ( i\j-r^'&r 
K, 



- (-l)"sin(V2^Kc<^c)cxp (30) 
3b) the extended singlet 



(-1)" cos(x/27ri^,^e) exp ( i^j —d, ] (31) 



"2^ 



3c) and the triplet pairing 



sin(V27r^(^s) exp ( iJ — i?, 



(32) 



— (c^ ± 



V2 
exp 



cos (fi^i?. 
Sin {jfd. 



(33) 



Note that the smooth part in Eq. (|30|1 corresponds to 
the usual BCS-type pairing while the oscillating terms in 
P0|l and (|31|l describe the eta-pairing superconductivity 



2.2 Phases 

With the results of the previous section for the excita- 
tion spectrum and the behavior of the corresponding fields 
Eqs. p7|) ~ (|19|l we now analyze the weak-coupling ground 
state phase diagram of the model |^ (sec Fig.^. 

Let us first consider the case J7 = 0, where the weak- 
coupling analysis shows existence of two different phases: 
in the case of antiferromagnetic exchange, at Jxy > 0, 
there is a gap in the charge excitation spectrum while 
the spin sector is gapless. Ordering of the field (pc with 
vacuum expectation value {ipc) = leads to a suppression 
of the CDW and superconducting correlations. The SDW 
and Peierls correlations show a power-law decay at large 
distances [221 • Due to the C/(l)-spin symmetry, K* > 1 
and the "in-plane" correlations dominate in the ground 
state, 



+ (-l)V-i/^'' (34) 



while the longitudinal spin correlations 

{S'[r)S'{0)) ~ + (-1)^-^= (35) 

and Peierls correlations decay faster. 

We now focus on the case of ferromagnetic exchange 
between spins. Ki U — Q and J^y < there is a gap in 
the spin excitation spectrum while the charge excitation 
spectrum is gapless. As common in the half-filled band 
case, the gapless charge excitation spectrum opens a pos- 
sibility for the realization of a superconducting instability 
in the system. Moreover, due to the J7(l)-symmetry of 
the system, ordering of ips with vacuum expectation value 
{(fis) = ^/tt/8Ks leads to a suppression of the CDW and 
singlet correlations as well as 5~ = ±1 channels of the 
triplet pairing. However, the 5*^ — 5*^ and triplet correla- 
tions in the S'^ — channel show an identical power-law 
decay 

(5^(r)5^(0)) = {0!,sir)OTsm (-1)^-^ (36) 

at large distances and are the dominating instabilities in 
the system. 

Let us now consider the effect of the on-site Coulomb 
repulsion. At J^y > the easy-plane antiferromagnetic 
phase remains unchanged at [/ > 0. However, at J^y < 
the TS+SDW^^-' phase is stable only towards influence of a 
weak U < — Jxy on-site coupling. In the case of repulsive 
Hubbard interaction, at [/ > — Jxy a charge gap opens. 
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This regime corresponds to the appearance of a long-range 
ordered antiferromagnetic (Neel) phase 



(S'^(r)5^(0)) - constant 
in the ground state. 



(37) 



Xy-Phase 

Ac ^0 
As =0 



(LRO) 

Ac ^0 
Af, 




TS + SDW'-^') 



/t 



Ac J^O 
As =0 



Fig. 1. The weak-coupling phase diagram of the itinerant 
XY-model at half-filling. Solid lines indicate borders between 
phases. The shaded region qualitatively marks the transition 



into the XY magnetic phase {J. 



3 Numerical results 

We use the density-matrix renormalization-group (DMRG) 
method |37I38| to study the ground-state properties of 
this model. Our calculations have been performed for open 
chains up to 120 sites using the infinite-size version of the 
DMRG routine. A comparison with the finite-size algo- 
rithm, which requires more CPU time and memory, does 
not give a substantial improvement of the results. For 
most of the numerical results reported here we have kept 
400 states in each block, which produces truncation errors 
smaller than 10"^. 

In order to reduce edge effects we average correlation 
functions C{\i — j\) over a number of pairs (i, j) of lattice 
sites separated by the same displacement r :— \i — j\ |^ . 
Typically we take nine pairs and for each value r we place 
the pairs as close to the center of the chain as possible. 

The asymptotic behavior of correlations (e.g. expo- 
nents) has been determined by an appropriate fitting of 
the data gOI. 



3.1 Excitation spectrum at [/ = 



2.3 The ferromagnetic transition 

Let us now discuss the ferromagnetic transition in the itin- 
erant XY model {U = 0). The very presence of this transi- 
tion can already been seen within the weak-coupling stud- 
ies, however detailed analysis of the phase diagram close to 
transition is out of scope of the continuum- limit approach. 
As we obtained, at Jxy < 0, \Jxy \ 'C t, the charge excita- 
tion spectrum is gapless and the spin excitation spectrum 
is massive. However, in the limit of strong ferromagnetic 
exchange \Jxy\ ^ t, the model is equivalent to the XY 
spin chain. Therefore, with increasing coupling one has 
to expect a transition from the regime with massive spin 
and massless charge excitation spectrum into a insulating 
magnetic phase with gapless spin excitations. On the other 
hand, in the case of antiferromagnetic exchange Jxy > 
the weak-coupling study shows a phase with gapless spin, 
gapped charge and dominating easy-plane spin correla- 
tions. One expects that this phase evolves smoothly to 
the strong coupling limit. 

The Jxy <-> —Jxy asymmetry is already seen on the 
level of the Hartree regularization of the band- width cut- 
off parameter W = 27rt as given by the Eqs. H14f) 



W. 



eff 



2n 1 



'Jxy 
2TTt 



(38) 



The weak-coupling approaches fail when the effective di- 
mensionless coupling constant \gi\ = = ~ 1. This 

condition immediately gives Jx'y = —irt. As we show be- 
low, using the DMRG studies of chains up to L = 120 
sites, indeed the transition into the ferromagnetic easy- 

-it. 



plane ordering discussed above takes place at J; 



xy 



Let us start from the limiting case of the itinerant XY 
model {U = 0) and analyze its excitation spectrum. The 
charge and spin gap for a half-filled L-site system are eval- 
uated by 



L 
^'2 



L 



L 



1 +Eo 



L 
^'2 



'L L 
2'2 



(39) 
(40) 



respectively, where Eq {N^ is the ground-state energy 
for 7V| up-spin and down-spin electrons. The extrapo- 
lation for L ^ oo is then performed by fitting a polyno- 
mial in 1/L to the calculated finite-chain results. Figure |21 
displays the extrapolated values as a function of Jxy 

We observe the following four sectors: at Jxy > the 
system is characterized by gapless spin and gapped charge 
excitation spectrum, while the weak-coupling ferromag- 
netic sector exhibits gapless charge and gapful spin de- 
grees of freedom. 

Moreover, our numerical results show the presence of 
two new regions. At j!fy^ « — 3t a charge gap opens, while 
the spin gap starts to decrease and finally closes at jiy'^'' sa 

< /"'^ 

xy ^ fJxy 



—At. This defines two new sectors: for Jilf ■* < J. 



both the spin and charge sectors are gapped, while at 
Jxy < Jxy^ the spin sectors become gapless. There are no 
indications for further transitions in the system. Note that 
similar behavior of the gaps, with interchange of spin and 
charge degrees of freedom, was first observed by Sikkema 
and Affleck in the Penson-Kolb model |4H . 
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Fig. 2. Spin and charge excitation spectrum of the itiner- 
ant XY model at half-filling. Depending on the behavior of 
the gaps four sectors can be distinguished. The approximate 



boundaries J^y^' 
lines. 



-3t and J, 



(c2) 



-At are indicated by grey 
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Fig. 3. The longitudinal (circle) and transverse (square and 
diamond) spin-spin correlations in the case of strong easy- 
plane ferromagnetic J^y = —8t (sector I) and antiferromag- 
netic Jxy = 8t (sector F) exchange. 



3.2 Correlation functions at U = Q 

To investigate the nature of ordering in the different phases 
we study the behavior of the correlation functions. In the 
sectors with gapless excitation spectrum and at half-filling 
we expect the usual expression for correlation functions 

C{r) = {O\r)O{0)) - Air-^' + (-l)M2r-''= (41) 

consisting of a smooth part decaying with exponent 9i 
and an oscillating part decaying with 62- In determining 
the asymptotics of correlation functions we focus on the 
dominating part given hy 9 = min{0i, 02}- 

In the following we will present results for correlation 
functions in different sectors of the phase diagram. 

3.2.1 Sectors I and F [Ac ^ 0, As = 0): The XF-phases 

In Fig. 121 we have plotted the longitudinal and transverse 
spin-spin correlations in the case of strong easy-plane ex- 
change. Although the amplitudes of the transverse cor- 
relation functions are different, the estimated exponents 
are similar. In the case of ferromagnetic exchange we ob- 
tained 6 w 0.57, whereas for the antiferromagnetic ex- 
change we have 9 « 0.61. The results are in a good agree- 
ment with the exact value 9 = 0.5 obtained for the XY- 
model 021 • The longitudinal correlation functions decay 
faster. The calculated exponents 9 « 1.79 (for Jxy = —8t) 
and 6 w 1.66 (for Jxy ~ 8t) are close to the exact XY- value 
9 = 2. 

The asymmetry of this model is clearly seen in Fig. 01 
where the ground state energy as a function of J^y is pre- 
sented. As we observe from Fig. 0] in the case of ferromag- 
netic exchange the ground state energy of the itinerant 
model becomes very close to that of the spin-i XY chain. 




Fig. 4. Ground-state energy Eo of the XY-model and the 
itinerant XY-model as a function of coupling Jxy/t for a half- 
filled L = 100 chain. 



3.2.2 Sector II {Ac = 0. Ag ^ 0): The TS+SDW^^) 
regime 

Let us now focus on the case of ferromagnetic exchange 
Jxy < at [/ — 0. The bosonization results predict a 
suppression of the CDW and singlet correlations, whereas 
SDW^^-* and triplet correlators show identical power-law 
decay (cf. with Eq. (|^ ). Furthermore, they are the domi- 
nating instabilities in this phase. Figure displays DMRG 
results for the singlet- and triplet-pair correlation func- 
tion. One can clearly observe a strong triplet-pair cor- 
relation. Note that the on-site and extended singiet-pair 
correlations show an almost identical behavior. This is ex- 
pected from the bosonization results H30() and (|31|l since 
the smooth part of the on-site singlet correlations ()30|l 
does not contribute due to In the double logarithmic 
plot (see lower figure) all correlation functions indicate a 
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tions exhibit an algebraic decay in which the transverse 
spin-spin correlation function decays faster. The calcu- 
lated exponent of the longitudinal spin correlation func- 
tion is, in agreement with bosonization results, close to 
that of the triplet-pairing correlations. Compared to the 
other results, the scaling behavior of the spin correlations 
is not very good in this sector. However, we have verified 
that with increasing system size and number of states kept 
in the DMRG algorithm the region with algebraic decay 
increases. Nevertheless, the numerical estimates for the ex- 
ponents are less reliable than those for other correlation 
functions. 

To complete the weak-coupling picture of sector II, we 
performed calculations for the density-density correlation. 
The results are shown in Fig. |3 Since in the double loga- 



10 



30 



Fig. 5. Pair correlation functions in the case of ferromagnetic 



exchange Jx 



-2t and [/ = (sector II). The lower figure 



shows the algebraic decay of the triplet and singlet correlation, 
plotted on a double logarithmic scale. 



power-law decay with fast decaying singlet-pairing corre- 
lators {9 2.18) and a slowly decaying triplet correlation 
function {6 « 1.03). The results are in a good agreement 
with those predicted by bosonization. 
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Fig. 6. DMRG results for the longitudinal and transversal spin 
correlation function, plotted against the real space distance 
\i — j\ (upper figure) at Jxy = —2t (sector II). The exponents 
(lower figure) were calculated using a suitable subset of the 
data to reduce finite size effects and numerical inaccuracies at 
large distances. 



In Fig. El we show calculations for the longitudinal 
and transversal spin-spin correlation for ferromagnetic ex- 
change {Jxy = — 2t). We observe that the correlation func- 
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Fig. 7. DMRG results for the density-density correlation func- 
tion at Jxy = —2t (sector II) including the average value of this 
correlation which removes the even-odd-r oscillations. 



rithmic plot we observe strong oscillations we additionally 
calculate the average value 031 

C(r) = i [C(r - 1) + 2C{r) + C{r + 1)] (42) 

to smoothcn the curve. As its clearly seen from the lower 
part of Fig. the oscillations are removed, but the es- 
timated exponent remains almost unchanged. Thus the 
DMRG result indicates a fast decay of density-density cor- 
relations, in agreement with the bosonization results. 

Therefore, we can conclude that coexisting triplet-pairing 
and antiferromagnetic SDW*-^^ ordering are the dominat- 
ing instabilities in this sector. 



3.2.3 Sector III {Ac ^ 0, As ^ 0): The intermediate phase 

In this subsection we analyze the intermediate sector at 
—At < Jxy < —it which is absent in the weak-coupling 
phase diagram (cf. with Fig. 
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Wc start to examine the asymptotic behavior of the 
superconducting and spin-spin correlations. In Fig. |Hlwe 
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tudinal spin correlations show well-established long-range 
order. 

The appearance of LRO is consistent with the U{1) (E) 
Z2 spin-symmetry of the present model The continu- 
ous U{1) symmetry is generated by the operators and 

, while the discrete Z2 symmetry comes from the invari- 
ance with respect to the — > —S^ transformation. Since 
the SDW^ ordering violates the discrete Z2 and transla- 
tion symmetries, the true LRO state is not forbidden. 



-0.01 

10° 




10 20 30 40 



Fig. 8. Pair correlation functions at Jxy = —3.5t and U = 
(sector III). The lower part shows a double logarithmic plot 
with linear fits. 



3.3 Behavior for nonvanishing on-site interaction 

3.3.1 Excitation spectrum at U ^ 

Let us now consider the effects of a repulsive Coulomb in- 
teraction on the ground state phase diagram of the model. 
We start with the excitation spectrum. 

From the bosonization results we know the general ef- 
fect of the Coulomb repulsion on the phase diagram which 
displays itself in an enlargement of the charge gap sectors 
at the expense of the spin gap sector. Fig. I lUI shows charge 
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Fig. 9. Spin correlation functions at J:cy = —3.5t and U — 
(sector III). The lower part shows a double logarithmic plot 
with a linear fit. 

present DMRG data for the pairing correlation functions. 
As is clearly seen from the figure, especially from the log- 
arithmic plot, the superconducting correlations decay ex- 
ponentially in agreement with the presence of a charge 
gap. 

In Fig. iniwe plot the spin-spin correlation functions. 
The logarithmic plot shows that the transverse spin corre- 
lation functions decay exponentially. Contrary, the longi- 
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Fig. 10. Spin and charge excitation gaps of the itinerant XY 
model ai U = t (upper figure) and U = 2t (lower figure). 

and spin gaps for U = t and U — 2t. One can clearly see 
that sectors I and F, where we have a finite charge gap 
Ac > 0, are enlarged. As a consequence the spin-gapped 
phase (sector II) becomes smaller with increasing U and 
finally vanishes completely. Already at U = 2t the charge 
gap is always finite. Thus the main effect of the presence of 
Coulomb interactions is the suppression of sector II, i.e. a 
reduction of the region with dominating superconducting 
correlations. In analogy with the U — case we conclude 
that the sectors with magnetic correlations become dom- 
inating. 
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3.3.2 Correlation functions at [/ 7^ 



In the following we analyze the effect of the Coulomb in- 
teractions on pair and spin correlation functions. We will 
focus on the behavior in sectors II and III where Ac = 0, 
As ^ and Ac 7^ 0, As 7^ 0, respectively. 

In the TS + SDW^^^ phase we consider the coupling 



J xy 



-1.5t Sit U ~ t as a. representative point. The phase 



is characterized by a spin gap of magnitude As ~ 0.1 3t 
and massless charge mode. The asymptotic behavior of the 
pair and spin correlation functions is plotted in Fig. II Hand 
Fig- El respectively. One can clearly see that the triplct- 
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Fig. 11. Pair correlation functions in the ferromagnetic phase 
at Jxy = —1.5t and U = t (sector II). The lower figure shows 
the algebraic decay of the triplet and singlet correlation, plot- 
ted on a double logarithmic scale. 



pairing and longitudinal spin-spin correlations represent 
the dominating instabilities in the system. Unfortunately 
the accuracy of the numerics is not sufficient in this case 
to verify that the exponents are still exactly identical. In- 
stead, we find 6 ~ 1.2 (triplet pairing) and 6 w 1.39 (lon- 
gitudinal spin). 

In the SDW^^' phase with LRO we compute the cor- 
relation functions at U = 2t and J^y = —2t. The presence 
of a charge gap Ac ~ 0.38i leads now to an exponential 
decay of superconducting correlations. 

On the other hand, as is clearly seen from Fig. ^1 the 
longitudinal spin-spin correlations show a true LRO while 
the transverse spin correlations decay exponentially. 

In ferromagnetic phase we use as a representative point 
U = 2t and Jxy = —At. As one can observe from Fig. [TH 
the longitudinal spin-spin correlations are exponentially 
suppressed, while the decay of the transverse ferromag- 
netic spin correlations is almost identical to that of the 
standard Xy-chain. 
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Fig. 12. Spin correlation functions at Jxy ~ —1.5t and U = t 
(sector II). Lower figure is plotted on a double logarithmic scale 
showing an algebraic decay for the longitudinal correlator and 
exponentially decaying behavior for the transversal correlation 
function. 
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Fig. 13. Spin correlation functions at Jx 



-2t and U = 2t 



(sector III). Upper figure indicates LRO fiuctuation for the 
longitudinal spin correlation function. 



4 Conclusions 

Motivated by recent experimental findings that show evi- 
dence for the competition or even coexistence of supercon- 
ductivity and magnetism we have investigated the ground 
state properties of an itinerant XY model. 

First we considered the case of vanishing on-site Coulomb 
interactions. The behavior of spin and charge gaps as func- 
tion of the spin-coupling J^y allows to distinguish four 
different phases (cf. with Fig. I15() . For antiferromagnetic 
interactions Jxy > the spin gap vanishes, but the charge 
gap is always finite. The observed behavior of the correla- 
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Fig. 14. Spin correlation functions at Jxy = — 4t and U — 2t 
(sector I). 



tion functions indicates a smooth evolution to the limiting 
case of spin-1/2 antiferromagnetic XY chain at J^y oo. 

For ferromagnetic couplings Jxy < three different 
phases appear. Already for weak interactions a spin gap 
opens, but the charge sector is gapless. Here SDW'^^ and 
triplet correlations, which decay with similar power-laws, 
are dominating, i.e. this regime exhibits a coexistence of 
antiferromagnetic ordering and triplet superconductivity. 

At Jif^ w -3t the spin gap is maximal and a charge 
gap opens. This intermediate phase, that extends up to 

(c2) 

Jxy ~ — 4i, shows long-range order in the longitudinal 
spin correlation, whereas superconducting correlations are 
suppressed and decay exponentially as expected for the 
case of a finite charge gap. Finally, at Jxy > — 4t via a 
spin-flop transition the system again enters a XY phase 
characterized by vanishing spin but finite charge gap. Here 
the behavior is similar to the ferromagnetic XY model. 



The presence of a repulsive on-site Coulomb interac- 
tion U has a strong effect on the phase diagram. Generi- 
cally it leads to an enlargement of the sectors with non- 
vanishing charge gap at the expense of the sectors with 
spin gap. Already aXU = 2t the charge gap is finite for all 
values of the exchange coupling Jxy Therefore the phase 
where antiferromagnetism and triplet superconductivity 
coexist is no longer observed and magnetic correlations 
become dominant everywhere. Only for small values of 
the Coulomb interaction there is still a finite window of 
coexistence possible. 

Fig. ^] summarizes our findings. The phase diagram 
shown combines results from bosonization (cf . Fig. ^ and 
DMRG. However, it is difficult to determine the location 
of the phase boundaries numerically. E.g. it still has to 
be clarified whether the coexistence phase extends up to 
the XF-phase such that two LRO SDW*^^-* phases exist 
(Fig. [T3i). Alternatively, only one LRO SDW^^^ phase 



LRO SDW^^^ 
# 




Xy-Phase 



/ LRO 




Fig. 15. The possible ground state phase diagram of the itin- 
erant XY-raoAe\ at half-filling. Solid lines mark second or- 
der phase transitions between the phases as obtained from 
bosonization. The dashed line corresponds to the spin-flop 
transition from the LRO SDW*-^-* into the ferromagnetic XY 
phase [Jiy^^ ~ — 4t). The dashed-dotted line marks the metal- 
insulator transition from the spin-gapped metallic phase with 
identical decay of triplet superconducting and SDW'^' correla- 
tions into the LRO antiferromagnetic (SDW*-^-*) phase {Jx'y'^ ~ 
—3t). Our numerical results do not exclude the possibility 
that the multicritical point given in 



A 



is an artifact of the 
bosonization approach. In a more realistic scenario of the phase 
diagram the SDW^^' always separates the superconducting 
phase from the easy-plane ferromagnetic phase as it is shown 



B 



exists such that it always separates the superconducting 
from the ferromagnetic XY-phase (Fig. llSbL 
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